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Abstract 

A new non-perturbative approach to quantum theory in curved 
spacetime and to quantum gravity, based on a generalisation of the 
Wigner equation, is proposed. Our definition for a Wigner equation 
differs from what have otherwise been proposed, and does not imply 
any approximations. It is a completely exact equation, fully equivalent 
to the Heisenberg equations of motion. The approach makes differ¬ 
ent approximation schemes possible, e.g. it is possible to perform a 
systematic calculation of the quantum effects order by order. An iter¬ 
ative scheme for this is also proposed. The method is illustrated with 
some simple examples and applications. A calculation of the trace of 
the renormalised energy-momentum tensor is done, and the conformal 
anomaly is thereby related to non-conservation of a current in d = 2 
dimensions and a relationship between a vector and an axial-vector 
current in d = 4 dimensions. 

The corresponding “hydrodynamic equations” governing the evolu¬ 
tion of macroscopic quantities are derived by taking appropriate mo¬ 
ments. The emphasis is put on the spin-| case, but it is shown how 
to extend to arbitrary spins. Gravity is treated first in the Palatini 
formalism, which is not very tractable, and then more successfully in 
the Ashtekar formalism, where the constraints lead to infinite order 
differential equations for the Wigner functions. 
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Introduction 


One of the greatest problems we are facing in physics today is gravitation: 
as is well-known we do not possess a well-behaved theory of quantum grav¬ 
ity. Even semi-classical gravity, i.e. quantum held theory in curved space- 
time, where the gravitational held is kept as a classical background, is not 
that well-behaved. It seems that perturbation theory is just not applicable 
to problems involving gravitation, hence new, non-perturbative approaches 
have to be found. In this paper, inspired by recent results in QCD (especially 
the study of quark-gluon plasmas), we propose a non-perturbative approach 
based on the Wigner operator. The equation governing the behaviour of this 
operator is completely equivalent to the Heisenberg equations of motions for 
the helds. This Wigner equation is the quantum analogue of the classical 
Boltzmann or Vlasov equations, and thus allows for a better intuitive 

understanding of the problems. 

Start with a hat spacetime manifold and a Dirac spinor held, on it. 

Instead of just looking at 'ijj{x) or its Fourier-transform 'ip{k) we consider a 
“distribution” W{x,p) which is a functional of '0 and depends explicitly on 
both X and p. This “distribution”, the Wigner operator, is, however, a non¬ 
positive operator, which is why we put the term “distribution” in quotation 
marks. The explicit form is given in the next section; it is essentially a kind 
of Fourier transform of the density matrix. The Dirac equation for -0 gives 
rise to an integro-diherential equation for W, which is readily generalised 
to curved spacetimes. This allows us to attack the problem of fermions in 
curved space-time from a new angle, and introduce new, non-perturbative, 
approximation schemes. 

The next step is naturally to write down a similar equation governing the 
gravitational helds. Since fermions causes torsion p|, |^, it would at hrst 
sight seem convenient to introduce two Wigner operators and for 

the spinor connection, and the vier-bein, e“, respectively. These Wigner 
functions will not however be covariant, thus introducing an unwanted gauge- 
dependency. In the Ashtekar formalism instead, we can introduce two co¬ 
variant Wigner functions, one for the S'[/(2)-connection, and one for 
its conjugate Ef. In both cases, however, do the equations become highly 
complicated, and we wont be able to say much about them. We will basi¬ 
cally, in this paper, consider the background geometry as hxed (semi-classical 
quantum gravity). Let us note that such an approach should hnd many ap¬ 
plications: Hawking radiation, the early universe etc. We will not, however. 
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deal so much with these applications in this paper. 

This is not the hrst time the powerful Wigner-equation methods have been 
suggested as a framework for quantum helds in curved spacetime. Some 
earlier work by Calzetta et ah and Kandrup, 0, uses either Riemann nor¬ 
mal coordinates (and hence is only applicable in a sufficiently small region 
of spacetime) or a mode-decomposition (which only works for static space- 
times), they do, however, show the connection between this kinetic approach 
and the usual path-integral or Green’s function-based approach; the Wigner 
function can be expressed as a Fourier transform of a two point function with 
respect to the middle point between the two points (the precise relationship 
is given in the next section). In general, neither the Fourier transform nor 
the concept of a middle point makes sense in a general curved background, at 
least not in a coordinate independent manner. Therefore one could attempt 
to use Riemann normal coordinates in order to be able to dehne x — x' for 
X, x' points on the manifold, or one could make use of a mode decomposi¬ 
tion of the solutions to the equations of motion of a free held to generalize 
the Fourier transform. Neither of these two approaches are completely gen¬ 
eral, although they do represent useful calculational short-cuts. Winter, 0, 
has proposed an approach similar to the one put forward here in which one 
integrates along curves (geodesics, to be precise). But his approach needs 
the solution of as well the geodesic as the geodesic deviation equations, a 
daunting task in general. This again essentially restricts the usefulness of his 
method to Riemann normal coordinates or to other approximation schemes 
(hrst quantum corrections to the classical kinetic equation). The solution of 
the geodesic equation is not needed for the approach presented here, even 
though we at an intermediate stage uses parallel transport along geodesics. 
Halliwell, and GellMann and Hartle has shown that the Wigner function 
appears naturally in a quantum theory of histories, [Q, as for instance quan¬ 
tum cosmology; it is the only other place in which the Wigner function of the 
gravitational degrees of freedom is studied. It should hnally be mentioned 
that Fonarev too, has extended the Wigner function technique from hat space 
time to a curved background, |]^. Like the approach proposed here, he uses 
the structure of phase-space as a cotangent bundle, our result, however, is 
much closer related to the analogy with Yang-Mills theory. He furthermore 
emphasises the case of the scalar held, whereas this paper mostly studies 
spinor helds, with some comments on the gravitational held. The derivation 
of the conformal anomaly and of the “hydrodynamic equations” are also new. 
Most of the calculations are done in d = 4 dimensions for concreteness but 
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are valid for all d. In some particular cases the case d = 2 is also considered. 


The Method 


We start out by considering the spacetime geometry as being hxed, i.e. we 
begin by studying QFT in curved spacetime. Our main source of inspiration 
is the papers by Vasak, Gyulassy and Elze |]^, in which the Wigner operator 
formalism is derived for Dirac helds interacting with Yang-Mills helds. We 
will give a brief introduction to their results here, as this will indicate how 
to generalise to curved spacetime. 

Consider a second quantised Dirac held D(a;), the Wigner operator is dehned 
to be, [|, I 

W{x,p) = j e + -y) ® 'ipix - 


note that W takes values in the Clifford algebra, a point we’ll use later. This 
dehnition can be made gauge invariant by replacing 

= e-^^y^i){x) ( 2 ) 


by the gauge covariant generalisation 

i’ix- ]^y) = 


(3) 


where V is the (gauge) covariant derivative. Writing = a: ± we can 
write 


W{x,p) 




d^y 

Wr 


(4) 


where 


U{b, a) = Fexp 




A^{z)dz 



(5) 


where denotes the Yang-Mills held, and the path of integration is a 
straight line: z{s) = a -|- (6 — a)s, which goes from a to 6 as s goes from 
0 to 1. The Dirac equation then implies the following equation for W: 


[Pi, + 


W = -igYX.W 


( 6 ) 
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where is an integral operator involving the held strength tensor: 


where we have introduced the triangle operator 


1 - 


s 


giisA 





(7) 



( 8 ) 


It is understood that the derivative with respect to momentum always acts 
on the Wigner distribution alone and not on any of the other terms. 

One should also mention the relationship to Green’s functions. If G{x, x') de¬ 
notes the Green’s function, G{x, x') = ('0(a;)'0(a;')), then the Wigner function 
can be rewritten as 

IG(a;,p) = (IG(x,p)) oc j e~^^'PG{x - ^y,x + ^y)d‘^y 


i.e. as a Fourier transform of the two point function with respect to the 
distance between the two points. This shows the relationship with the more 
familiar approach to quantum theory using Green’s functions. 

The curvature 2-form appears through the holonomy, i.e. through re¬ 
expressing the derivative of the parallel transporter in the principal bundle 
U{b,a) in terms of a integration along a closed curve. The above equation 
actually holds for being an operator, but for a classical background one 
can pull the held strength tensor out of the integral over y to obtain 

where joGi are the usual spherical Bessel functions, jo Ji/ 2 )ii ~ Jz/ 2 - We 
will now analyse this equation and the steps leading to it. 

First the obvious translations: the gauge held is replaced by the spinor 
connection = oj^aahi whereby the gauge covariant derivative becomes 
the spinor covariant derivative, and the held strength tensor becomes the 
Riemann-Ghristohel curvature tensor Let us next consider the path of 
integration. In hat space we integrate along the straight line from = x—^y 
to = a; -|- i.e. we go in the direction ±y from x for a “period of time” 
|||?/||. The curved space-time analogue is now obvious: y is a tangent vector. 
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and we integrate along the (nniqne) geodesic with tangent ±.y at x, the length 
of the curve segment at each side of x is ^||y||. Note that this simplihes things 
a lot: a priori we would have expected y to he in the manifold, but now we 
see that it really lies in the tangent space. Hence for each given point x in 
the manifold we integrate over all possible tangents y lying in the tangent 
space at x, T^M, which is a flat space. This means that we keep the measure 
■4-^. To summarise: in flat space-time x,y,p all belong to the same space, 
whereas in curved space-time x is in the manifold, y in the tangent space 
and p in the cotangent space. The dot-product y ■ p is then just the pairing 
between T^M and T*M. What we have done is to make extensive use of the 
intimate relationship between Yang-Mills theory as a theory of connections 
on a principal bundle and general relativity as related to connections on the 
frame bundle, [ 0 , |, 

We will follow customary notation and denote flat space indices by Roman 
letters from the beginning of the alphabet, and curved spacetime ones by 
Greek letters. In curved spacetime, the Clifford algebra relation { 7 ^, 7 ^,} = 
implies that the Dirac matrices are in general x-dependent, we intro¬ 
duce the vierbein as the transformation connecting them to the usual flat 
(constant) Dirac matrices 0, |[: 



7m = e“7a (10) 

The connection between “curved” and “flat” indices is established by this 
vierbein. The covariant derivative then reads, |§], 

= + ( 11 ) 

with (Tab = |*[ 7 a, 76 ] the generators of the Lorentz algebra so(3,1) in the 
spinor representation, and is the analogue of the generator of the gauge 
algebra in the Yang-Mills case. The spin connection 0 ;“^ is related to the 
Christoffel symbol through 





V Ce 


C p p,!^ 


ca^bQ 


7] e 




( 12 ) 


and it is this quantity which is the analogue of the Yang-Mills potential 
Even though we deal with fermions, we will in general assume that no torsion 
is present, i.e. Tjjj, = T^^. This is a valid assumption since torsion do not 
propagate and since we only want to deal with helds in a given classical 
background. 
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With these comments the Wigner operator equation now reads 


m + Y\ e^Pa + 


W = --K^X^W 


( 13 ) 


where now 




d 


with the triangle operator given by 


A = _V 

^ a ^ IJ- 

OPa 


(14) 


(16) 


The curly brackets denote an anticommutator. Both the curvature two- 
form and the Wigner operator are Clifford algebra-valued, thus the Clifford 
algebra and not only so(3,1), appears here as the analogue of the internal 
symmetry group of Yang-Mills theory. One should remember that so(3,1) ~ 
spin{3, 1) C C(3,1), so we have in a sense “extended” the gauge algebra from 
the Lorentz algebra to the full Clifford algebra, with the price, of course, that 
we then no-longer have a Lie algebra but a Clifford algebra instead. Also 
note that so(3,1) ~ spin{3, 1) acts naturally on the Clifford algebra C(3,1) 
through X H->■ [x, cr“^], it is under this Lie algebra that the Wigner function 
transforms. 

Here we have just made the obvious translations from Yang-Mills theory to 
minimally coupled Dirac fermions in a curved background. Non-minimally 
coupled fermions, i.e. including a coupling to torsion, should be dealt with 
by making the substitution where is the 

torsion (the antisymmetric part of the Christoffel symbol). This is because 
the curvature only enters through the commutator [V^^, Vy] = i?“^crab (via the 
holonomy) and this gets modihed to the above mentioned linear combination 
of curvature and torsion when the latter is present, [§. In a similar fashion 
one could introduce gauge degrees of freedom, coupling to a Yang-Mills held 
Aj by adding igAp:'^ to V^, where T*, are the generators of the gauged Lie 
algebra. We would then have to add igFji^Tk to the right hand side of the 
commutator [V^,Vi/]. The resulting Wigner function and its equation of 
motion would then also be gauge covariant. 

Covariance is ensured by noting that -0 transforms as 0 ^ Uip where U is the 
transformation matrix {U = exp{iaab{x)a°'^) for the purely gravitational case 


7 





and U = exp(zQ;ab(a;)cr“^ + iga^{x)Tk) with a coupling to a Yang-Mills field), 
and that transforms covariantly (adjoint representation) US/ 

one then sees W —^ UWU~^. 

The Wigner equation is then an infinite-order differential equation. It is this 
equation which is the subject of this study. 


Clifford Decomposition 

The Wigner operator takes values in the Clifford algebra, since it is given as 
a product of two spinors. Hence it can be decomposed. In d = 4 a basis for 
the Clifford algebra is given by 1, 75, 7 a, 7 a75, as is well-known. We thus 

write 

W = 5 + P 75 + V“ 7 a + ^“7a75 + T“Vab (16) 

where then, of course, S,V are a scalar and pseudo-scalar, vector 

and axial-vector and an antisymmetric tensor. A similar splitting-up of 
the operators can be made using the properties of the Dirac matrices. The 
left hand side can be written 


+ 7 '^ (^Pd + 

= [m + + l5YBd] W 




W 


(17) 


having used 


^dcba^ 


7 cr = i{r]°''^r] - ) 7 a + 4e 757 a 

which follows from the definition of aab and the standard trace formulas for 


the Dirac matrices, see e.g. Itzykson and Zuber []W 
We will write the right hand side as 




(18) 


where then contains all the curvature information and no Clifford algebra 
information (i.e. it is proportional to 1). 

With this notation we get 


mP + AaM“ - HaV“ = -VacVbdr^^A'^ 
mV^ - A^S + B^V - i{r]efS^g - r]eg6j)A^Tfs 







mA'^ + A^V - B'^S + i(77e/<^3 - VegSj)B^rf3 

= -4e\J-^^S + 8z{vaAt-VacSt)r^^V 

mr! + - 

where we have used the symmetry properties of i.e. to 

remove some terms involving rjA. 

One will often be able to assume V = A°' = 0, in which case the equations 
can be simplihed a little bit. 

The equations involve not only the curvature 2-form but also its dual, since 
eabcdJ''^’' = -^me^(jo - iji)-^{*Rf^)Vbd (19) 

with 

* ■= 2^abcdRfy ( 20 ) 

One knows that in Yang-Mills theory self-dual and anti self-dual curvature 
2 -forms play an important role (they correspond to instantons, iJ, 0,0]), 
it is therefore encouraging that the dual of the curvature 2-form also appears 
in this gravitational setting. 

One should also note that the combination r]abSfJ°‘^^ is proportional to the 
Ricci tensor, rjabS^j^ = ^7r~^RK-^Rjf{jo — iji), which implies that this 
combination vanishes for vacuum solution of general relativity, i.e. in Ricci- 
flat spacetimes = R^f'iJd^'^ = 0. This implies that in many backgrounds 
of real physical interest, the right hand sides of these coupled equations sim¬ 
plify. 

Another important consequence of this Clifford algebra decomposition con¬ 
cerns the classical limit. In this limit one would expect the Wigner function. 


P_ ) =51 + V%„ 


considered as a 4 x 4 matrix to be diagonal: 


/ A 


w = 


A 


9i 


\ 


92 / 


F+ 

0 
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This leads to the following set of coupled equations 


mS - AoV° 

= 0 

( 21 ) 

5oV° 

= 0 

( 22 ) 

mV° - 

= tPabSlJ^^^S 

(23) 

B'^S 


(24) 


where, as we have seen, the right hand side in the third of these vanish if the 
spacetime is Ricci-flat, i.e. a solution to the vacuum Einstein equations. 

The second of these of equations imply Vo = 0, which then, from the hrst, 
yields m = 0 , hence we are left with the following pair of equations 

= -iVab^lr^’^S (25) 

B^S = (26) 


Thus in a Ricci-flat spacetime, the only quantum corrections enter through 
the second of these, as the hrst become AqiS = 0 , i.e. the quantum corrections 
will depend only on the dual of the curvature two-form. The hrst equation 
implies, then, that S is an eigenstate of the diherential operator + 
with eigenvalue — 2 ipQ, whereas the second equation determines the depen¬ 
dency on the other components of the momentum. We will return to the 
case W = S m. later sections. 

Let us note in passing that the extra term —2ipQ added to the Dirac operator 
can be interpreted as a coupling to an Abelian connection, or as correspond¬ 
ing to a Dirac operator in 5 dimensions (as in the families index theorem. 
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For later use, we also give the Clihord decomposition in d = 2 dimensions. 
Here 7 “ corresponds to <Ti,a 2 (the Pauli matrices) and both 75 and (Tab cor¬ 
responds to CTs = —iaia 2 = — <^ 2 ]. The covariant derivative can be 

written as = 5^, since there is no spinor connection in d = 2, [|^. With 
this (Latin indices going from 1 to 2) the Wigner equation reads 


m - a\pi + -idi) 


w = -jZKaXJRaf, i era, £-W [> (27) 


With the decomposition 


W = Wo 1 + Wi + W 3 


(28) 
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we then arrive at 


mWo - {pi + 

mWi - eijW^ - {pi + ]4di)Wo 
mVhs — i€ij{p^ + 


(29) 

-Uk 

(30) 

(31) 


We will return to this set of equations in the section concerning the trace of 
the energy momentum tensor in h = 2 and d = 4. For now let us just note 
that the set of coupled equations in d = 2 can be re-expressed as an equation 
for IFj alone by using the hrst and third to express IFo, IF 3 as some operators 
acting on IF*. Explicitly 


fRo = 

m 'idWj 1 , 

e^e^RapJ ^ + {Pi+^idi)W 

2 m opp m 2 

(32) 

fRs = 

in -dWj 1 , ,• l.ov'XTTPi 

pRiRapJ ^ ^ + eij{p + td)W^ 

4m opp m 2 

(33) 


leading to 

11 1 ■ 
{w?pij - meij)W^ - —{pi + -idi){pj + -idj)W^ 


1 


+ iriK€^^{pi + iridi) ( e'^R^pJ 


dpp ) 


\^r]’^R^elR^pR^sJ‘^ 


d'^W^ 


dppdps 


- - e“ e 7 Rap Rys 


dppdps 

1 1 '^dW^ 

--melejRapipk + 

(34) 


Remembering J = jo(|A) — zji(|A) and that A is hermitian, the hermitian 
and anti-hermitian parts of this equation must be satished separately, which 
then leads to two coupled equations which IR has to satisfy. 
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Example: Covariantly Constant Curvature 

When the curvature tensor is covariantly constant, Va-R^^ = 0^ can take 
the square of the equation for W 

(^m - jabc, = 0 


to obtain 

0 = - 275/1 ■ B]W + 

^2 a 2 

j!K|757.B“.e7''^/^{'^e/.M'}| (35) 

We can split this into two equations by taking the hermitian and anti- 
hermitian parts out. For the hermitian part we then get the constraint 
equation generalising the mass-shell condition 

^ K"^ 

0 = {a,„ w}} 

(36) 

whereas the anti-hermitian part becomes the kinetic equation proper (as in 
the Yang-Mills case, |I|) 

0= | 7 «. 4 “ + 757 ,B“.e; 7 ‘‘ 7 {Y{^^^,,, 7 }j (37) 


The mass-shell condition, (^), can also be seen as a momentum diffusion 
equation of the Fokker-Planck type ^ , but the important thing to note is 
that it is not on its own a kinetic equation specifying the state of the matter 
helds. Comparing this with a standard classical kinetic equation for a phase 
space distribution function JF, [Q 


d 

dpu 


(a^JF) 


1 

2 dp^dpu 




f{p,x)R 


(38) 
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we see that Oi,, the dynamical friction in onr case is 

a,^ = 0 (39) 

while the momentnm diffusion coefficient, is (simply take W = 5 to get 
the non-spin terms) 

K, = - %iK)R%Rt + spin terms (40) 

and f{p, x), a kind of momentum force term, dually reads 

f{p, x) = — (m^ — + Clifford algebra terms (41) 

and is thus an operator (remember: A p + d + u). 

Thus in this particularly simple case, we see that the mass-shell condition 
can be interpreted as a diffusion equation in momentum space. Now, this 
momentum diffusion is provided by the quantum and curvature corrections 
to the mass, i.e. the difference between and it is thus a pure quantum 
phenomenon. This is a kinetic interpretation of the uncertainty principle and 
of the renormalisation of the mass. It is intuitively pleasing to see that the 
momentum diffusion coefdcient is given simply by the curvature, and is thus 
a purely geometric object, even though its source is purely quantum. For a 
curvature which is not covariantly constant, the mass-shell condition contains 
higher derivatives with respect to the momentum, such terms do not have a 
direct classical interpretation and therefore represent pure quantum effects. 
The main difference between (pTl) and the result for Yang-Mills theory is the 
appearance of the anticommutator. As we have seen earlier, this anticom¬ 
mutator can be traced back to W being Clifford algebra-valued and that the 
generators of the Lorentz algebra for Dirac fermions is precisely Oab which is 
an element of the Clifford algebra. Remember that the elements ^alb of any 
Clifford algebra C{r,s) generate the Lie algebra spin{r,s) which is isomor¬ 
phic to so(r, s). For r = 3,s = 1 this is the Lorentz algebra in d = 3 -|- 1 
dimensions, while it for r = s = 1 is the corresponding algebra in d = 1 -|- 1 
dimensions, 0]. It is this interrelationship between the Clifford algebra and 
the Lorentz algebra (which is then the gauged algebra for gravitational sys¬ 
tems) which accounts for the anticommutator. 

In the extreme, classical limit W = S, the kinetic equation reduces to 
0 = 2i{r]afPde - r]aePdf)A'^e^‘^Wj^^S - 

OPiy 

( 42 ) 
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This can be seen as an equation involving only dynamical friction and not 
momentum diffusion. Writing it on the classical form 

(43) 

OPu 

where the “force term”, /i, contains the operator A, we see that the dynam¬ 
ical friction is given by 

= (riaiVde - Vd^m)(p‘ - - 

(44) 

It is very interesting to note that the resulting equations for W splits into 
two set of equations, one containing only momentum diffusion and the other 
only dynamical friction. In both cases are the sources for the processes given 
by quantum effects, i.e. vanish in the classical limit if one ignores spin effects. 
Furthermore, the coefficients a^, are in both cases given by the curvature 
and are thus geometrical quantities as one would expect intuitively. 

Again, in the general case we would get higher derivatives of momentum, 
the “diffusion” equation containing all even powers of and the “friction” 
equation all the odd powers. These extra terms are pure quantum effects 
and have no direct classical interpretation. 

To round off the discussion we also give the splitting of the Wigner equation 
in d = 2. As we have seen, in this case we can make do with an equation 
for Wi alone. The hermitian part of this equation is seen to be in the case of 
covariantly constant curvature (corresponding to taking J = 1) 


- meij)W^ 

1 If. 

-{PiPj + ^i{Pidj + pA) - = - g etelR^pAs 

+ ]^KekARAP^ + 


for the hermitian part, and 




(45) 

(46) 
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for the anti-hermitian part. Once again, we can interpret the hermitian part 
as a momentum diffusion equation, both this time we also have a dynamical 
friction contribution 

1 rP' 8 

The diffusion coefficient is 

><% = (47) 

while the dynamical friction becomes 

+ ]^idi) (48) 

which is then a differential operator in this case. The force term is again re¬ 
lated to the mass-shell constraint, although in a somewhat more complicated 
form 

- pV - - P^d^) + (49) 

The anti-hermit-an part is once more a pure dynamical friction equation 
(it would have a momentum diffusion term too, if the curvature wasn’t co- 
variantly constant), but there is no source term. Thus there is quite a big 
difference between the results in d = 2 and d = 4 dimensions, a difference 
which cannot simply be referred to the different Lorentz character of the 

“distribution” in the two cases (in d = 2 we considered a vector, and in 

d = 4 a scalar), but is very much due to the difference in the structure of the 
two Clifford algebras and the fact that there is no spin connection in d = 2. 
This is an important caveat. 

The A-Expansion: Quantum Corrections 

The case of covariantly constant curvature corresponds, as we have seen, to a 
kind of classical limit involving only few quantum corrections. In this section 
we will commence a more systematic study of quantum corrections. This is 
done by noting that A appears multiplied by h, and that Planck’s constant 
only enters in this combination. An expansion in A is therefore related to an 
expansion in h. Thus, one can calculate quantum corrections by expanding 
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the spherical Bessel functions in X to a given order. To do this we need the 
standard formulae 


io(^) 

ii(^) 


oo 

sin z . 


z 

sin 2 ; 


n=0 
COS 2; 


■ir 


v2n 


(2n + 1)! 




OO 

5:( 

n=l 


-ly 


2 n 


^2n-l 


(2n + 1)! 


Now, A is an hermitian operator and the Bessel functions only appear in the 
combination Jq — iji, consequently it is rather simple to make a separation 
into hermitian and anti-hermit-an parts of the Wigner equation. One of these 
will then only contain odd powers of A and the other only even powers. 

Due to the complicated nature of the operators appearing in the Wigner 
equation, we cannot simply take the square and compare with the classical 
kinetic equations. Instead we have to make do with the equations on the 
“Dirac form” (in contrast to the “Klein-Gordon form” resulting from taking 
squares), which do not have a direct classical analogue. The interpretation 
will therefore not be as precise as one would perhaps have wanted. 

Writing 

OO 

W{x, p) = J2 {x, p) (50) 

n=0 

an remembering = |h.A, we can collect terms with the same power of 
Planck’s constant. Doing this we get 

m + e((7“(p^ + ^fV^) 

for the “classical” contribution, hP, which is recognised as the same as for 
the case of covariantly constant curvature treated in the previous section, 
whereas the first correction satisfies 

m + e^ 7 “(p^ + 



( 52 ) 


= 


2 ^ ‘ 


be 


A. 

dpt 


■W^^\aA + 


= --we^7“i?; 


>bc 

fiiy 


A. 

dp^ 


TP®, 


a 


be 


(51) 
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while the second order contribution satisfies (remembering that the vierbein 
is covariantly constant) 


In general we will have a recursive scheme 

» + + iv,.) + imeST-fljL <7^,1 

= terms involving only with k < n (54) 

The terms on the right hand side will involve more and more momentum 
derivatives of the Wigner functions and similarly higher and higher or¬ 
der covariant derivatives of the vierbeins and curvature two-forms. Higher 
order derivatives of curvature can in general be considered as related to fluc¬ 
tuations of the geometry, thus the higher order Wigner functions are de¬ 
termined by the fluctuations of spacetime, moreover they couple to higher 
order derivatives with respect to the momentum of the lower order Wigner 
functions, terms which do therefore not have a classical interpretation (only 
first and second order momentum derivatives appear in classical kinetic equa¬ 
tions). This shows how the expansion in A is closely related to pure quantum 
effects with no classical analogue. One cannot, however, guarantee that this 
expansion is equivalent to the standard loop expansion in quantum field the¬ 
ory. In general there will be no such simple relationship. 

One should also take notice of the fact that the momentum derivative op¬ 
erator is symmetric in its indices, this implies that only the symmetric part 
of Vp...Vo- acting on the curvature will contribute. From the commutator 
relation we get 

V(^V.) = - R^laab (55) 

whence it follows that the right hand side of the recursive scheme not only 
includes higher and higher order derivatives of the curvature but also higher 
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and higher powers of it. 

If we introduce the operator Y by 

VW- m + + + W (56) 

then we can write the recursive scheme as 

yH,'(") = F<"> (57) 

where — 0 and F*"* depends on k < n. Thus, on the formal level 

U’*">(i,p) = V'“"F<">(x.p) , n>l (58) 



(59) 

(60) 
(61) 

(62) 

( 63 ) 
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for the scalar part, 


mV + - B-Va + 

2 ^ opu 

for the pseudo-scalar contribution, 

- A^S + B^V - - PegSj)A^Tf3 + 

Ki/B'-m - [-<& - 26V “V + 

for the vector part, while the axial vector contribution turns out to be 
mA'^ + A^V - B'^S + - peg5j)B^Tf3 _ 

“A-T'o - '-iKe'irf^A [- 4 e‘^^„S + Si(Stii - S^S^V] 

and hnally 

- 

[i26V 

for the tensor part. On the right hand side of the recursion relation we have 
terms of the form 



Now, here the covariant derivative is to act on the two form i.e. 

^ r>bc o r)bc , pA rybc , pA pbc 

'' P^gu ^P^gu ' ^ pg-^Xu ' ^ pu^Xg 

and hence Vp doesn’t involve any Clifford algebra elements {a^c is only the 
generator of so(3,1) in the spin | representation). Therefore we can move 
the (constant) Dirac matrix 7 “ outside the covariant derivation. A Clifford 
decomposition of this term is thus straightforward, and we obtain 

2 d'^A^ 
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p 

V 




he 




^‘"dp.dpp 

1 (9^ 9 


^^'dp^dpp 3 


K6V 


d^V 

dpudpp 


+ 


rjbfVcee^i'^pR] 


be ^ ^ q-ef 


A 


3 -, 

2 (9^ 9 

- 3 ev:(vxya^ + 


6c 




d'^V 

dpudpp 


2 Ke>^€\. ^ ^ - 


''a^ be ^ dg 


>^^'dp,dpp 


To have a look at a solution we can take the extreme case 
VpW(x,j9) = So{x,p) = 

As we have seen earlier, this is only possible provided m = 0, in which 
Sq has to satisfy the coupled set of equations 

- = 0 
-B^So + A5„ = 0 

The second of these yields upon insertion of the explicit form for Sq 

u;l^ = -KR%{a'^ppp + n 

which implies = 0. Inserting this into the hrst we arrive at 

Pd + ^'i(^d{-^tJ-(^''^Pi'Pp + dp^'^Pv + dpi) = 0 
which is only possible if 

= 0 

(3'" = 2iA' + const. 

7 = const. 


(64) 

case 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 
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Thus 5o is of the so-called Jiittner form [Q 


exp{/3{x){fM{x) -ifU^{x)) 


where (3 is the inverse temperature, /r the Gibbs energy and the four 
velocity, i.e. (3^ = f3U^, 7 = f3fi. Let us now make a similar Ansatz for the 
hrst quantum correction, i.e. fT^) = ^ 5 ^ = A/”'exp(—-f j3^p^ + 7 ). 
The equations then reduce to 


d 


1 (9^ 




flUJ 


dpudp 


(72) 
-5o (73) 


then evaluating the differentiations with respect to the momenta and collect¬ 
ing powers of these, we arrive at the following set of conditions 


Pd 


d^a'^P = 0 
\ieP{dt^^lP^ + \i^ePRlia'^Ppp = 0 


= 0 

^eP{ul^6tVcd + ^^R:t)S^ = ^zKeP{V,R^J,)aP^So 
from the first of the equations for iSi, while the second give us 


(74) 

(75) 

(76) 

(77) 


Combining these we get the following set of conditions 


(78) 


d^oAP 

RlWP 


0 

0 


2 pd 


In a general spacetime the hrst two of these will give = 0 as for the 
lowest order term. Hence the solution iSi will once more be of the Jiittner 
form, albeit with a much more complicated expression for P{x)Up{x) = [3^. 
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To get more information about the recursive scheme presented in this section, 
we can hnd the hermitian and anti-hermitian parts of the squared equations 
(i.e. the Wigner equation on “Klein-Gordon form”), as this is where we 
expect to see the analogy with classical kinetic theory most clearly. 

We essentially calculated the square of Y in the previous section, the square 
of is similarly 


(^(l))2 


— 

36 


d2jY(0) 

dpudpp 



dpndpe 



( 79 ) 


A Clifford decomposition of this is complicated by the 7 “ part of and the 
quadratic appearance of We note that (Fd))^ ig hermitian, hence it 

only contributes to the momentum diffusion equation generalising the mass- 
shell constraint, whereas the kinetic equation proper (which comes from the 
anti-hermitian part of F^) does not get any contribution from lower order 
terms. The momentum diffusion equation then reads 


(m^ -A^-B^-2A- F75)iyd) 


16 “ P^dp^dpo 

ra2pp(o) 
t ^Pu^Pp 




36 




I 

\ dp^dp, ’ 

( 80 ) 


At the next order, O(h^), however, there will be contributions from Wd) 
to the kinetic equation proper for Vhd) too, namely 



a^wd) 

dpudpp 



dp„dp^,dp^ 



In general, the momentum diffusion equation will contain only squares, whereas 
the kinetic equation will contain only cross products (always in the form of 
an anticommutator) of the lower order Wigner functions. 


The Conformal Anomaly in ci = 2 and d = A 

Let us calculate the trace of {T^^), denoted by (T). The non-vanishing of 
this quantity is the conformal anomaly when m = 0 , [@. 
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The general expression for (Tab) in terms of the Wigner function is given by 

{Tab{x)) = TT [ -faPb{W)d^p (81) 

Jt*m 

It follows that (in d dimensions, W = {W)) 


(r):=(W) = Tr [ rPaWd'^p 

Jt*m 


Now, from the equation for W 


(82) 


^ + e;(V^) 


W = XW 


it follows that 


7 '^PaW = - m)W - XW 


Hence 


(T) = Tr 


'T*M L 


- m)W - XW 


d!^p 


We also know 


d 


XW = ——iW X (other terms)) 
OPu 


(83) 

(84) 

(85) 


i.e. that XW is a total derivative, and thus that its integral over the cotan¬ 
gent space at x vanish, leaving us with 


(T) = Tr [ {Urey^-m)Wd^p 
Jt*m 2 


( 86 ) 


From this we see that (T) measures the failure of W to satisfy a Dirac 
equation with mass 2m (thus it is most interesting when m = 0). This 
equation is valid for all d. To carry out the trace we need the Clifford 
decomposition, which is then d dependent. 

For d = 2 we get 


{T{x)) = l-i[ ey’’df,Wbd‘^p-m [ W^p (87) 

2 Jt*m Jt*m 

where we have written W = Wq 1 -|- Wa cr“ - 1 - IF 3 with a, 6 = 1, 2. In¬ 
troducing the current density {ja{x)) := f Wad^p and the number density 
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(n(x)) := f WqcPp we can write this in terms of these macroscopic qnantities 
only as 

{T{x)) = - m (n(x)) ( 88 ) 

For massless helds, the conformal anomaly is then related to the non-conservation 
of the cnrrent ja, snch that (T) 7 ^ 0 if and only if ja is not conserved in ex¬ 
pectation valne. 

In d = 4 we similarly get 


(r(x)) 





- - m{n) 


ct^p 

(89) 


where we have introdnced the vector cnrrent {ja) '■= f Vad^P, the axial cnr¬ 
rent (ka) := / Aad'^p and the nnmber density (n) := J Sd'^p. We see that in 
d ^ 2 the possible existence of a conformal anomaly or not is not as simply 
related to the qnestion of the conservation of a cnrrent as in d = 2 . 

The fact that an anomaly can be expressed in terms of the Wigner fnnction 
integrated over the cotangent space, snggests a closer relationship between 
W and anomalies, in particnlar with the spin complex, |1l2|, |T^]. The very 
natnre of the Wigner fnnction, or rather the entire Wigner-Weyl-Moyal for¬ 
malism where operators are replaced by symbols on the cotangent bnndle, 
makes the translation of analytical properties into geometrical or algebraic 
ones possible. Snch a translation is at the heart of index theorems, and it has 
in fact been shown that the Atiyah-Singer index theorem can be related to 
the classical limit of the Wigner-Weyl-Moyal (WWM) formalism for the ordi¬ 
nary Heisenberg algebra, |]^. The way in which this generalised WWM formalism 
relates to index theorems is presently under study. 


The Hydrodynamic Equations: Moments 

The previous section calculated the trace of the energy-momentum tensor from the 
Wigner equation, this is only one macroscopic quantity which one can define. Let 
us define the following moments of the Wigner function (for simplicity we suppress 
the (•)) 


n(x) 


Tr 



/ 


S(Ap 
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U^{x) 

:= Tr 

1 p‘^W(T^p = J p°'Sd!^p 

ja{x) 

;= Tr 

l‘xaWd‘^p= J Vad'^P 

ka{x) 

;= Tr 

j IblaWd^p = J Aad^p 

Tab{x) 

;= Tr 

1 PalbWd‘^p = J PaVbd'^p 

Sab{x) 

;= Tr 

^ aabWd'^p = J Tabd'^p 

Tab{x) 

;= Tr 

PalblbWd‘^p = J PaAbd*p 

X{x) 

:= Tr 

XbW d^p = J Vd^p 

Xa{x) 

;= Tr 

1 PalbWd^p = J PaPd^p 

^abci.x') 

:= Tr 

^ PaCTbcWd'^p = j PaTbcd^p 


These are all the zeroth and first moments of the Wigner function in d = 4 di¬ 
mensions. Most of these quantities have a direct physical interpretation, n{x) is 
the number or energy density, is a momentum density, ja a current, ka and 
axial current, Tab the energy-momentum tensor and Sab the spin density, while Tab 
is a kind of “pseudo-energy-momentum tensor”, x and Xa are related to chirality 
and \abc represents a spin-momentum interaction term. The equations of motion 
for these macroscopic quantities constitute the corresponding set of (quantum) 
hydrodynamic equations, which are derived by simply taken zeroth and first mo¬ 
ments of the Wigner equation in its Clifford decomposed form. Hence each of the 
five equations in this Clifford decomposition gives rise to two equations for these 
moments. If we write 

Aa=Pa + \iDa Ja^^ = (90) 

we get the following ten equations 

(91) 

(92) 

(93) 


mn-T - -^iDaf + Bak°- = 0 
mUb - J PbPaV^^cfp + ^WaTb'' + BaTb'" = e\, 
mx + T+ ^iDak'" - Baf = 0 
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mxb +J PaPbA^d^p+^iDaTi,°‘ - ( 94 ) 

m/ - [/'^ - ^ib‘^n + B'^x 

+ + = 0 ( 95 ) 

- JPhP’^Scbp - ^ib'^Uh + B^^Xh 

-iiVefSg - VegSj){ j PhP^'T^U'^p+]^ib'^ 

+Ae,j^B^X^^ = -i6tJe^,R%eln 

+2e\/je^,R%elx 

-^il‘^‘'r]bfVceJe^R%els^^ ( 96 ) 

mk‘^ + X'^ + - B^n 

-Ae,f^{X^f^ + \tb^s^b = 0 ( 97 ) 

mrjf + I PhP'^Vd^P + \ib\h - B^^Uh 

+i{r,,f5i - Veg6j)B^X/^ - 46 ,^/( Jp.p^Rf^d^p + Ub^X,!^ = ^e\Je^,R%eln 

-Si{5t5i - 5-5i)Je^,R%elx ( 98 ) 
ms^f _ = 0 ( 99 ) 


+ \ilPhP^^V^^dS - Ib^^T^/ 


^bKJ^ = 

2 1^1 

+ i^c^''^db ~ 


sy>,^)je'iit^y^k‘‘ (100) 


where T = T^,t = r“. Some of these have direct physical interpretations, e.g., 
as we saw in the previous section (^) gives a kinetic expression for the confor¬ 
mal anomaly when m = 0. In analogy with this, we will refer to the equation 
for r, (1^), as the pseudo-conformal anomaly for lack of a better word. Another 
important equation is ( p5[ ) which states that the momentum density and the mass 
times the current density (which is also the velocity density) are not identical, the 
difference between the two is an indication of heat flow. Also, equation (^) gives 
an expression for the spin density in terms of the velocity/current density ja and 
the axial current ka- 

The quantities without a direct physical interpretation such as Tab,X^Xa, ^abc can 
be eliminated from these equations. 
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One should also note that we are not able to eliminate the second moments appear¬ 
ing in this set of equations. In general one will get an infinite hierarchy of moment 
equations. This can be truncated by brute force at any given stage resulting in a 
set of approximate hydrodynamic equations. The second moments have a physical 
interpretation in terms of viscous pressure, since 


J php'^W(fp = J i 


:= Tr / php'^Wcfp = I php'^Sd^p 


( 101 ) 


is the viscous pressure tensor, [^. The remaining second moments are then higher 
Clifford algebra analogues of this, lacking a straightforward classical interpretation. 


Arbitrary Spins 

Clearly the fermions are the most difficult to treat due to their Grassmannian 
nature and that is why we have chosen to consider such fields in great detail. 
Arbitrary spins are not that difficult. Consider a field equation for a field <I> of 
some spin s, 

{D, + = gj 

where Dg is some differential operator (first order whenever s is half integral and 
second order when s is an integer), Mg is a mass-term including couplings to 
curvature such as for s = 0 and for s = 1, s' is a coupling constant and J 
is a source term, as we have seen above such an equation gives rise to an equation 
for the associated Wigner function by the “minimal substitution” 



and the “source term transformation” 

(103) 

where X is some integro-differential operator containing the source. For bosons 
with spin s, A is a 4s tensor and IF is a 2s tensor, for s = 1/2 we have already 
seen that X carries a Clifford algebra index (i.e. two spinor indices) and three 
Lorentz indices, two of which must be understood as coming from the Lie algebra 
so{r,d — r), for s = 3/2 we would get a Wigner function which had two more 
Lorentz indices, and X would also have four extra Lorentz indices to account for 
the vector index on the Rarita-Schwinger field. It is clear that for s > 1 the 
notation quickly becomes cumbersome. 

All Wigner functions are maps from the cotangent bundle into some algebra, for 
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a scalar field the target space is simply C, for spin-i, it is C{r, d — r), for s = 1 we 
get {T*M ( 8 ) 0 )( 8 ) (r*M ( 8 > 0 ^) with g the internal gauge algebra and so on. In general 
fermions will have Wigner functions with values in C{r, while 

bosonic fields will have Wigner functions taking their values in gauge 

degrees of freedom are handled by simply enlarging the target space by tensoring 
it with 0 ( 8 ) 0 ^. 

A spin-0 field, 4>, would then give rise to a Wigner equation of the form 

+m^ + iR 

while the Maxwell field gives rise to 

(105) 

and a Rarita-Schwinger field would give rise to 

(e>e + IR.p = (106) 

where we have omitted the A and ip subscripts on the Wigner functions for the 
cases of s = 1, 3/2 and suppresses the spinor indices in the latter case. 

Yang-Mills fields could be treated analogously but will have more complicated ki¬ 
netic equations due to the non-linearity of their field equations. 

Internal degrees of freedom is also treated easily. The fields will now be cross 
sections in some associated bundle, and the parallel transporter has to include not 
just the curvature effect but also the connection in the corresponding principal 
bundle (the gauge field). One simply replaces the covariant derivative with the 
appropriate gauge covariant derivative containing the gauge field and the metric 
connection. If the field transforms in the representation p of 0 then the Wigner 
function will transform in the representation p 0 p, A p is the fundamental repre¬ 
sentation, then this is the adjoint one. 

The Gravitational Field: Palatini Formalism 

Next we write down the equations for the gravitational field. This can be viewed 
from two points; ( 1 ) as a means of calculating the back-reaction of the quantum 
fields on the space-time geometry (keeping the energy-momentum tensor fixed), or 
( 2 ) as constituting full-fledged quantum gravity (giving a set of coupled equations). 
We immediately face a problem: all-though the connection is the formal analogue 
of the Yang-Mills fields, we cannot take over the results by Else, Gyulassy and 


IT^ = y(o)w^ (104) 



Vasak, as the Einstein equations are first order in the connection, while the Yang- 
Mills equations are second order. In the presence of torsion the field equations for 
gravity can be written 

elRf, = K(r; + ie“r) (107) 

S^ab = + (108) 

where T“ is the energy-momentum tensor, T its trace, the torsion and is 
given by 

Cab = ^iE'4)'y^aabil^ (109) 

with E = det e“, see Ramond Q. We will have to introduce two Wigner operators, 
one for the vierbein and one for the spin-connection: 

rjr = (110) 

l;1 ^ ( 111 ) 

Unfortunately these are non-covariant, as the spin connection transforms in an 
affine way under local so(3, l)-transformations. This non-covariance will then in¬ 
troduce an unwanted dependency on coordinate choices (gauge dependency). 
Finding the Wigner equations for these is straight-forward (albeit tedious), and 
will not be done here, they will have the form 

+ ( 112 ) 

“ (M ^ Z^) = “ (/^ ^ ^) (113) 

^ ( 7 ^ (J) 

Where Xy ’ ^ are integral operators containing the source terms. Expressing these 
sources in terms of the Wigner operator for the matter fields, e.g. 

Tab =< TV( J ^aPbWd^p) > (114) 

would then lead to a set of coupled integro-differential equations, constituting the 

full set of equations for quantum gravity. It is to be remembered that for a 
quantum field, we cannot take the curvature two-form out of the integral in the 
original expression for X. Therefore, the equation for the matter Wigner function 
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becomes much more complicated, in Q the analogous situation in QCD is treated. 
We will not attempt this level of generality here, just note that the straightforward 
generalisation from Yang-Mills fields to gravitational fields does not give tractable 
kinetic equations, at least not in the Palatini formalism. 


The Gravitational Field: Ashtekar Variables 


As Elze, Gyulassy and Vasak, Q, have developed the Wigner function technique 
for Yang-Mills fields, one would suspect that an approach similar to theirs can be 
fruitful if one uses the Ashtekar formulation of gravity. In this formulation, |15], 
one has a complex 5f7(2)-connection, A = Afdx’'o'a, where Ua are the generators 
of SU 2 and i = 1,2,3 is a spatial index. This formalism is based on the fact 
that the Lorentz algebra so(3,1) in four dimensions (and four dimension only) 
is isomorphic to the complexification of SU 2 - The Ashtekar formalism therefore 
makes use of two vital aspects of general relativity: the correct dimensionality and 
the correct signature of spacetime. Canonically conjugate to the connection, we 
have the “electric field”, E, 


i.e. 




( 115 ) 


Introducing the field strength tensor PT-, the constraints can be written as 


A = 

H = 

Ga = V,E^ = 0 


In exact analogy with the Yang-Mills case as presented in Q, we could define 


^ab 


{x,p) = 

Jt^m 


d y -ip-y / -ly-V pa 

{2nrf ' 


f: 




with Tq- = Ef. Considering the Hamiltonian nature of the system, it is, how¬ 
ever, more appropriate to introduce two slightly different Wigner functions (i.e. 
essentially splitting the above candidate into two), namely 


■pa6 

ijkl 




f ^ y r-^p y 

Jt^t, (2vr)3 

f ^ y r-^p y 

Jt^t. (2vr)3 








(116) 

(117) 


where we have used the global hyperbolicity of spacetime always assumed in a 
Hamiltonian formulation of gravity, M ~ S x M, and furthermore made explicit 
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use of the complex nature of the connection and the “electric field” (whence the 
daggers on the right hand side). 

The constraints can then be used to derive relationships between these two Wigner 
functions. The diffeomorphism constraint Di can be rewritten 

(118) 

Noting that the Wigner functions precisely the Weyl transforms of 

these two tensor products, we see that the quantum version of this constraint reads 

= Q (119) 

where * denotes the twisted product (the non-commutative product induced on 
the set of phase-space fnnctions by the non-Abelian product of operators), 

f * g ^ = fg+^if A g... (120) 


with f A g := {f,g} being the Poisson bracket. 
Similarly the Hamiltonian constraint H becomes 


nj _ ^ab M 

ri = e 


r,f-L/r 


, * VI, * TT', 

T bb' ij'i'j' 


= 0 


( 121 ) 


One should note that while the classical constraints Di , H are purely algebraic 
in the quantum versions become infinite order differential equations in 

the corresponding Wigner functions ~ equivalently, they become ^-algebraic, i.e. 
deformed. The last constraint, the Gauss one, however, gives rise to a differential 
eqnation for Lfj, analogous to the equation for the Dirac-Wigner function (or 
more precisely, to the Yang-Mills Wigner function as derived by Elze, Gyulassy 
and Vasak, citeEGV). 


Group Theoretical Arguments 

We can round off this discussion with a few comments about the algebraic structure 
behind the entire WWM-formalism. The flat space Wigner function has a natural 
group theoretical interpretation, |^. Define the operator 

n(u, v) = exp{iup — ivq) (122) 

then n(ri, v) form a ray representation of the (Abelian) group of translations in 
phase-space, 

n(n, v)U{u',v') = + u\v + v') (123) 
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This operator then gives the Weyl transformation, mapping an operator into a 
function on the classical phase-space 


Aw{u,v) = Tr n(ti,n)A 


(124) 


The Wigner function is the (symplectic) Fourier transform of the Weyl transform 
of the projection operator (for a pure state) 


i.e. 


W{x,p) = J ™^^lJ\Il{u,v)\'^|J)dudv 


(125) 


What we have done in this paper is to replace the Abelian group of translation on 
a flat phase-space with the non-Abelian group of parallel transport on the curved 
phase-space T*M. This splits up into two parts, the parallel transport in the base 
manifold M, which is generated by the momentum operator, and the (Abelian) 
translations in the fibre T*M, which is generated by q. 


In a previous work, [ 14 1, 1 have shown that one can generalise 11 to a very large class 
of algebras, most notably finite dimensional Lie algebras, their corresponding loop 
and Kac-Moody algebras as well as super Lie algebras and C^-algebras. To put it 
in algebraic language, then, what one does when going from flat space to curved 
space is to replace the usual Heisenberg algebra by the curved space analogue 


\PtJ.,Pu] = -Riiv 9 "'] = q''] = 0 

where R^u is the curvature two-form (not the Ricci tensor). 

The twisted product is given by 

Aw * Bw ■= {AB)w = Tr n(u, v)AB 


(126) 


(127) 


and can be written in terms of a kernel as |11, 


{f* 9 ){u,v)= / K{u,v,u',v',u",v")f{u',v')g{u",v")du'dv'du'dv" (128) 


with 

K{u, V, u', v', u", n") := Tr n(u, v)U{u', n')n(u", n") (129) 

Thus the quantity H is the essential ingredient in any generalised “Wigner-Weyl- 
Moyal formalism” or “symbol calculus”, Hence what is needed in general is 
(i) a “symbol map” sending operators (typically pseudo-differential operators) into 
functions on the cotangent bundle, and (ii) a map connecting two fibres, this latter 
map is simply the symbol of the parallel translator. In principle one could have 
a different map in (ii), but the symbol of the parallel translator is the simplest 
choice. 
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The formula we have derived for W was based on the phase-space of a clas¬ 
sical mechanical system, namely T*M, the Wigner function took its values in 
C{3, 1) ®T{T*M) then, where C(3,1) is the Clifford algebra and r(r*M) denotes 
the set of cross sections of the cotangent bundle. We used the group of parallel 
transport and translations along the fibres to generalise the Wigner function from 
quantum mechanics in flat spacetime. One can proceed to quantum field theory 
by means of second quantisation, although this is not usually easy to define in a 
curved background, [|[. Locally this always makes sense, but there are often global 
obstructions. Another method, which is not very common, is to treat the field and 
its conjugate momenta as the fundamental phase-space (which is then infinite di¬ 
mensional), the Wigner function is then a functional of these fields. This is very 
similar to the way one treats quantum gravity in quantum cosmology, where one 
then considers the wave function of the universe. The Wigner function has already 
been extended to this situation. 

Let us consider a field theory in Hamiltonian formalism and denote the fields and 
their conjugate momenta by 0, tt respectively. These can be either bosonic or 
fermionic. We then want a functional W such that 

(A) = J W[4>,7r]A[4>,TT]T>(j)T>TT (130) 

which generalises 


(A) = Tr pA = J Aw{x,p)W{x,p)dxdp 


If we have free fields, then the phase space is flat, and the parallel transporter 
becomes simply 

J uij>-vndx (131) 

From the quantum mechanical relation 

W{u,v) = Tr pll{u,v) = {Il{u,v)) 

which is the symplectic Fourier transform of the Wigner function, we conclude 

1 


W{u,v) := (n(u, u)) = 


e*^n(u, v)V(j)V7r 


Z{0) 


(132) 


Noticing the concrete form of H we see that exp(f5)n is equivalent to adding a 
source term (with the sources denoted by u,v respectively), hence 


W{u,v) = 


Z{u, v) 


(133) 
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The proper Wigner functional is then a symplectic (functional) Fourier transform 
of this quantity. This makes it possible to interpret the partition function (with 
sources) as essentially the Wigner functional of the vacuum state (or more general, 
the vacuum to vacuum transition). 


Conclusion 

We have generalised the work by Elze, Gyulassy and Vasak to gravitation and 
quantum fields in curved space-time. Thereby we obtained a set of exact equa¬ 
tions for QFT in curved space-time and even Quantum Gravity, which allows us 
to make non-perturbative calculations in these cases. The major draw-back is the 
complicated nature of the equations - especially in the full quantum gravity case. 
But on the other hand, we can develop a recursive scheme for these quantum 
corrections (the A-expansion), and we can therefore avoid perturbation theory 
all together. One should also note that this approach is on the one hand inti¬ 
mately related to the algebraic structure of the space of quantum observables (the 
canonical commutator relations for the fields), and on the other to the topological 
structure of spacetime (the Wigner function is, for Dirac fermions, a mapping from 
the cotangent bundle into the Glifford algebra, whereas for bosons it is a mapping 
from the cotangent bundle into a the tensor algebra). It was this interplay that 
allowed us to find dynamical expressions for the conformal anomaly in d = 2 and 
d = 4. 

We also saw that the kinetic equation satisfied by the Wigner function could be 
split up into two, one being the mass-shell constraint giving the quantum and cur¬ 
vature induced corrections to the mass, while the other was the kinetic equation 
proper. In certain classical-like situations these could be written as two Fokker- 
Planck equations, the one with no momentum diffusion the other with no dynam¬ 
ical friction. We derived expressions for these kinetic quantities. We also saw how 
quantum corrections modified this simple situation. 

By taking appropriate moments of the Wigner equation we arrived at a set of cou¬ 
pled equations governing macroscopic quantities such as energy-momentum tensor 
and current and spin densities. These were the corresponding hydrodynamic equa¬ 
tions. We saw that they, besides giving the kinetic interpretation of the conformal 
anomaly, also lead to an expression for the heat flow. 

Gravitational degrees of freedom was attempted handled first in the Palatini for¬ 
malism, in which we had to introduce Wigner functions for as well the vierbein as 
the spin connection. These was not, however, covariant, and the resulting equa¬ 
tions were too complicated. We then turned to the Ashtekar variables, where 
we could either introduce one Wigner function for in the standard way, or 
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use a canonical description to split this into two, one for F^- and the other for . 
These were by construction covariant. The constraints induced conditions on these 
Wigner functions, which were ^-algebraic, i.e. infinite order differential relations 
due to a quantum deformation of the product. 
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